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ABSTRACT 

We present kinematic parameters and absorption line strengths for three brightest 
cluster galaxies, NGC 6166, NGC 6173 and NGC 6086. We find that NGC 6166 has 
a velocity dispersion profile which rises beyond 20 arcsec from the nucleus, with a 
halo velocity dispersion in excess of 400 km s . All three galaxies show a positive 
and constant hi Hermite moment. The rising velocity dispersion profile in NGC 6166 
thus indicates an increasing mass-to-light ratio. Rotation is low in all three galaxies, 
and NGC 6173 and NGC 6086 show possible kinematically decoupled cores. All three 
galaxies have Mg2 gradients similar to those found in normal bright ellipticals, which 
are not steep enough to support simple dissipative collapse models, but these could 
be accompanied by dissipationless mergers which would tend to dilute the abundance 
gradients. The [Mg/Fe] ratios in NGC 6166 and NGC 6086 are higher than that in 
NGC 6173, and if NGC 6173 is typical of normal bright ellipticals, this suggests that 
cDs cannot form from late mergers of normal galaxies. 
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1 INTRODUCTION 

Brightest cluster galaxies (BCGs) are luminous ellipticals 
found near the centres of galaxy clusters and groups. They 
are further subdivided into giant ellipticals (gE), D galaxies 
which have shallower surface brightness profiles than gEs, 
and finally cD galaxies which have extended low surface 
brightness envelopes (Mathews, Morgan & Schmidt 1964; 
Schombert 1987, 1988). BCGs are often associated with 
powerful radio sources (Mathews et al. 1964), and extended 
cluster X-ray emission (Jones et al. 1979). In regular clus- 
ters, BCGs are invariably at the centre of the cluster po- 
tential, as defined by the cluster X-ray emission, galaxy 
distributions, and strong/weak gravitational lensing (e.g. 
Miralda-Escude 1995), and their velocities are usually close 
to the mean cluster velocity. BCGs are often observed to 
have extremely high globular cluster specific frequencies 
(e.g. Harris, Pritchet, & McClure 1995; Bridges et al. 1996). 

The origin of BCG halos is one of the most puzzling 
aspects of galaxy formation, and it is unclear whether they 
are formed at the same time as the galaxy proper, as pro- 
posed by Merritt (1984), or are accreted later. If accreted 



later they might form via dissipational processes, such as 
star formation in cooling flows (Fabian & Nulsen 1977), 
or by purely stellar dynamical processes, such as galactic 
cannibalism (Hausman & Ostriker 1978; McGlynn & Os- 
triker 1980), accretion of tidally stripped material (Rich- 
stone 1976) or "Galaxy Harassment" (Moore et al. 1996). 
For a review of BCG formation, see Tremaine (1990); see 
also Garijo, Athanassoula, & Garcia-Gomez (1997), and Du- 
binski (1998). 

The kinematics and metallicities of BCG halos provide 
important clues as to how they form. If they form in a dis- 
sipational process, then we might expect a low velocity dis- 
persion, as the energy of the material from which the stars 
form has been radiated away. If BCGs form in a dissipa- 
tive collapse, we would expect to see metallicity gradients 
in the halo (though metallicity gradients are not unique to 
dissipative formation models). If the halos form out of mate- 
rial from other cluster galaxies via dissipationless processes, 
then we might expect the halos to reflect the kinematics of 
the cluster as a whole, with high velocity dispersion and lit- 
tle net rotation; metallicity gradients would be shallower if 
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BCGs experience significantly more mergers than other clus- 
ter galaxies, since mergers are expected to weaken metallic- 
ity gradients. If the halos are primordial, and formed in a 
process involving violent relaxation, then the kinematics will 
reflect the underlying mass distribution. 

While there is considerable information on the central 
kinematics of ellipticals, there have been few studies of the 
halo kinematics at very large radii from optical spectroscopy. 
Most such studies have concentrated on normal ellipticals, 
and have found that the velocity dispersion profiles are ei- 
ther flat or decreasing with radius (e.g. Saglia et al. 1993; 
Carollo et al. 1995; Statler, Smecker-Hane, & Cecil 1996). 
Similar results have been found in the even fewer studies 
of BCGs (e.g. Carter et al. 1981, 1985; Tonry 1984, 1985; 
Heckman et al. 1985; Fisher, Illingworth, & Franx 1995). To 
date, one cD (IC 1101 in Abell 2029: Dressier 1979; Fisher et 
al. 1995; Sembach & Tonry 1996), and two dumbell galaxies 
(Abell 3266: Carter et al. 1985, and IC 2082: Carter et al. 
1981), have been found to have rising dispersion profiles. 

An alternative approach to the study of the kinemat- 
ics of the outer halos of cD and other galaxies is the use 
of globular clusters and planetary nebulae as tracers. Be- 
cause of the geometric constraints of current multi-slit and 
fibre fed spectrographs, this technique can be used only in 
the outer halos, and because of magnitude limits cannot 
be applied at distances beyond two or three times the dis- 
tance of the Virgo cluster. Nevertheless for M87 (Mould et 
al. 1990; Cohen & Ryzhov 1997) and NGC 1399 (Grillmair 
et al. 1994; Kissler-Patig et al. 1998; Minniti et al. 1998), 
globular cluster velocities have shown conclusively that the 
M/L ratio increases outwards in the halo, and thus that 
there are substantial dark matter halos in both galaxies. 
Arnaboldi et al. (1994) confirm this result for NGC 1399 
for a sample of planetary nebulae. Interestingly, Cohen & 
Ryzhov (1997) find that the globular cluster velocity disper- 
sion in M87 rises outwards, and in NGC 1399 the velocity 
dispersions of the planetary nebulae and globular clusters 
match smoothly onto the velocity dispersion of the Fornax 
cluster itself. These BCGs at least seem to merge smoothly 
into their host clusters, and the BCG and cluster dynamics 
seem closely linked (see discussion in Freeman 1997). The 
evidence from the distribution of lensed arcs in more dis- 
tant clusters points to the dark matter being more centrally 
concentrated than the X-ray gas (Miralda-Escude & Babul 
1995). In this case the dark matter will affect the dynamics 
of the cD halo significantly. 

We have embarked on a program to obtain very deep op- 
tical spectroscopy of BCGs. Here we present long-slit spec- 
tra for NGC 6166, NGC 6173, and NGC 6086, obtained at 
the 2.5m INT in La Palma in May 1996. The properties of 
these galaxies and their host clusters are listed in Table Jj| 
NGC 6166 is a luminous cD and a classic multiple-nucleus 
galaxy (though the 'nuclei' are probably not bound to NGC 
6166), centrally located in a rich cluster with a large cooling 
flow (~ 160 M /yr: Allen & Fabian 1997), and a regular, 
symmetrical X-ray appearance (Buote & Canizares 1996). 
Abell 2197 is 1.3 degrees north of Abell 2199 (Gregory & 
Thompson 1984), and is a sparser cluster with an irregular 
X-ray morphology dominated by two main concentrations 
(Muriel, Bohringer, & Voges 1996). NGC 6173 is a D galaxy 
found at the centre of one of these concentrations (subclus- 
ters), and appears to have a significant peculiar velocity (i.e. 



Table 2. Observing Log and Instrumental Setup 



Telescope 
Instrument 
Observing dates 
Camera 
Grating 
Dispersion 
Detector 
Resolution 
Spectral Coverage 
Slit Width 
Average Seeing 

Slit PA 
Exposure Time 



Isaac Newton Telescope (2.5 m) 
Intermediate Dispersion Spectrograph 
10-16/06/1996 
235mm 
R1200Y 
35 A/mm 
Tektronix 1024 CCD 
2 A 
4930-5730 A 
3.0 arcsec 
1.3 arcsec 
41/140/0 deg 1 
18.75/6.25/6.33 hrs 1 



1 For NGC 6166/6173/6086 respectively. 

relative to the cluster mean velocity). Abell 2162 is proba- 
bly best described as a poor cluster or compact group; we 
have not been able to find any X-ray data for this cluster. 
However, NGC 6086 is classified as a cD galaxy based on 
its surface brightness profile, and does not have a peculiar 
velocity. Thus, our three galaxies/clusters span a wide range 
of properties. 

NGC 6166 has previously been observed spectroscopi- 
cally by Tonry (1985) and Fisher et al. (1995) with 2 and 
3 position angles respectively (each include the major axis) . 
Unfortunately, both studies only extend to 10—20 arcsec 
from the galaxy centre, and over this limited range both 
find a flat or slightly decreasing velocity dispersion profile 
(we show below that this is consistent with our data, which 
extend much further out). Fisher et al. find that NGC 6166 
does not have significant rotation along either the major or 
minor axes out to 20 arcsec radius. We have only been able 
to find central velocity dispersions and line strengths for the 
other two galaxies (e.g. Burstein et al. 1987; McElroy 1995). 



2 OBSERVATIONS AND DATA REDUCTION 

Table ^ contains details of the observations. Spectra are in 
each case taken with the slit aligned with the major axis to 
maximise the signal to noise ratio. The spectra include the 
Mg and Fe absorption lines, as well as [OIII] and [NI] emis- 
sion lines (for NGC 6166 only). We also obtained template 
spectra of 28 stars, with spectral type ranging from G6 to 
Ml, and luminosity class III, IV and V. 

Figure |l| shows the central spectra for each of the three 
galaxies; for each spectrum, the best-fitting stellar template 
is also shown (see below), as well as the bandpass defini- 
tions for Mg2 , Fe527o , and Fes335 . Preliminary data reduction 
was done using standard FIGARO routines, and followed 
the procedures of Carter, Thomson and Hau (1998). This 
involves debiassing, flatfielding and wavelength calibrating 
the data. A slope in the spectrum due to detector rotation 
was corrected using the S-distortion correction routines in 
FIGARO. The template spectra were re-binned to zero ve- 
locity, using their published radial velocities, mainly from 
Wilson (1953). 

Following van der Marel & Franx (1993), the line-of- 
sight velocity distributions (LOSVDs) are modelled as Gaus- 
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Table 1. Properties of Target Galaxies and Host Clusters. We list data for our sample of three BCGs, and two others with rising velocity 
dispersions, IC 1101 in Abcll 2029 and the giant dumbell in Sersic 40/6 (Abell 3266). Column 1 designates the galaxy and its host 
cluster. Column 2 gives the galaxy morphological type from Schombcrt (1986), except for Sersic 40-6. Column 3 gives the total absolute 
blue magnitude assuming Ho=75 km s _1 Mpc _1 , and taking By values from Burstein et aJ. (1987), and Pclcticr et ai. (1990) for IC 
1101. Column 4 gives the central velocity dispersion: Burstein et aJ. (1987) for N6166, N6173, and N6086; Carter et ai. (1985) for Sersic 
40-6; Carter et aJ. (1981) for IC 2082; Fisher et ai. (1995) for IC 1101. Columns 5, 6 and 7 give the galaxy velocity, cluster velocity, 
and cluster velocity dispersion, respectively: Zabludoff et aJ. (1993a) for N6166/A2199 and N6173/A2197; Postman & Lauer (1995) and 
Zabludoff et aJ. (1993b) for N6086/A2162; Oegerle, Hill, & Fitchett (1995) for IC 1101/A2029; Teague, Carter, & Gray (1990) for Sersic 
40-6/A3266; Ellis et aJ. (1984) for IC 2082. Column 8 gives the Bautz-Morgan classification, and Column 9 gives the cluster richness. 
Column 10 contains information about the cluster's X-ray appearance: 'CF' means a cooling flow is present; 'Reg'/'Irr' means the cluster 
has a regular/irregular X-ray morphology. 
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Figure 1. The central galaxy spectra for NGC 6166, NGC 6086, and NGC 6173, in the rest frame of the galaxies. The lighter solid 
lines show the data, while the heavier solid lines are the best-fitting galaxy model for each galaxy. The scale is arbitrary. We have also 
indicated the definitions of the Lick Fe5270 and Fes335 bandpasses, and our redefined Mg2 index. 
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sians with small Hermite deviations, parameterised as mean 
velocity v, dispersion a, and hermite moments h n of order n. 
The parameters are recovered by the program kinematics 
written by H-W. Rix; for more details please see van der 
Marel & Franx (1993), Rix & White (1992), and Hau, 
Carter, & Balcells (1998). The Gauss-Hermite moments /13 
and hi are closely related to the skewness and kurtosis of 
the LOSVD respectively, and are therefore easy to visu- 
alise. Large /13 with opposite sign to the mean velocity is 
typical of rapidly-rotating systems, whilst for spherical two- 
integral models, /14 is proportional to the velocity dispersion 
anisotropy (van der Marel & Franx 1993). The higher mo- 
ments are more difficult to interpret. However, they are use- 
ful for identifying structures in the residuals of the model fit 
to the data. We measure and plot Gauss-Hermite parame- 
ters up to he,. 

The observed wavelength range is well known to give 
template mismatch problems (van der Marel & Franx 1993) . 
In order to minimize template mismatch, an optimal stellar 
template is employed for the direct fitting. This is estimated 
as the best linear combination of the 28 template spectra 
which when broadened with a Gaussian profile minimizes 
X 2 , the difference between the combination and the galaxy 
spectrum, after an initial guess of the velocity for each tem- 
plate and the dispersion, (Rix & White 1992). This optimal 
template is then used for recovering the LOSVD for the 
galaxy spectrum, modelled up to h&. This procedure is done 
for each of the individual spectra in each galaxy, thus reduc- 
ing radial template mismatch due to abundance gradients. 
The galaxy nucleus is located by a Gaussian fit to the light 
profile of the central few pixels. Monte Carlo experiments 
show that 2500 counts/pixel or higher gives satisfactory re- 
sults for the direct fitting (see Appendix A) . 

Our procedure allows the determination of an individual 
optimal composite template for each data point, To ensure 
that this procedure does not introduce systematic effects 
into our profiles we repeated the analysis for NGC 6166 us- 
ing the composite template for the central spectrum for all 
data points. There are no systematic differences in zeropoint 
or overall trend in any of the kinematic values derived. Sub- 
jectively, we estimate that the scatter is higher, justifying 
our decision to use individually determined optimal tem- 
plates. 

Emission lines could affect the measurement of 
absorption- line indices if they fall within one of the band- 
passes, and also can potentially introduce systematics in the 
absorption line profile analysis. Significant [OIII] and [NI] 
emission are observed in the nuclear region r < 11 arcsecof 
NGC 6166, but not in NGC 6173 or NGC 6086. Thus, for 
NGC 6166, pixels within a wavelength difference correspond- 
ing to 3<r from an emission line are excluded from the line- 
profile fitting. These regions are: 4988-5021A ([OIII]), 5179- 
5213A ([NI]), and 5397-5415A (strong sky line). Experimen- 
tation shows that excluding a few pixels has a negligible 
effect on the recovery of the LOSVD as kinematics are mea- 
sured over the entire spectrum. For spectra more than 11 
arcsec from the nucleus where no emission is present, only 
the region affected by the strong sky line is excluded, allow- 
ing fitting over the entire spectrum. Another difficulty with 
NGC 6166 is the presence of companion nucleus B, which 
is at a PA of 62°with respect to the central cD, and con- 
tributes significantly to the light in our slit between ~ 5—20 



arcsec from the galaxy centre along the NE side. Kinematics 
was not designed to recover two LOSVDs, so we restrict our 
analysis to the SW side of the galaxy. 

The parameters recovered by the Gauss-Hermite analy- 
sis are presented in Figures ^(a), (b), and (c), for NGC 6166, 
NGC 6173, and NGC 6086 respectively. 



3 STELLAR KINEMATICS 

3.1 Velocity Dispersions 

Central Velocity Dispersions 

We can compare our central velocity dispersions, 00 = 
320 ± 5, 295 ± 5, and 325 ± 8 km s -1 for NGC 6166, 6173, 
and 6086 respectively, with previous determinations. NGC 
6166 has been the best-studied of the three: Malumuth & 
Kirshner (1981), Tonry (1985), Burstein et al. (1987), and 
Fisher et al. (1995), find a = 350 ± 35, 303, 326, and 295 
± 6 km s -1 respectively. For NGC 6173, Burstein et al. 
(1987) and McElroy (1995) find <j = 261 and 240 km s -1 
respectively, both values lower than ours. For NGC 6086, 
Burstein et al. (1987) and McElroy (1995) find cr = 304 and 
325 km s -1 respectively. Overall, the agreement is good, 
and the present data are of much higher quality than any 
previous data. 

Velocity Dispersion Profiles 

From Figures |(b) and (c), we see that both NGC 6173 
and NGC 6086 have velocity dispersion profiles that decline 
outwards. Such behavior is found in most normal ellipticals. 
In contrast, Figure |^(a) shows that the velocity dispersion 
profile of NGC 6166 increases from ~ 325 km s -1 at the 
centre to ~ 450 km s -1 at 35 arcsec along the major axis. 
Tonry (1985) and Fisher et al. (1995) both have limited 
dispersion profiles for NGC 6166, and find a fiat and slightly 
decreasing dispersion profile respectively, out to ~ 20 arcsec 
along the major axis. This is consistent with what we find, 
since our dispersion only starts rising significantly beyond 
~ 20 arcsec. 

In Appendix A we analyse possible systematic errors 
introduced into the kinematic parameters by noise and im- 
perfect sky subtraction. We find that there is a bias towards 
low velocity dispersion when the signal-to-noise is poor, but 
this does not significantly affect the points plotted in Figure 
2(a-c). 

3.2 Rotation 

From Figure ^|(a), NGC 6166 shows modest major-axis ro- 
tation, and a weighted least-squares fit yields a slope of 1.07 
± 0.12 km s _1 arcsec -1 , or a rotation of ~ 45 km s -1 at 
40 arcsec. Tonry (1985) found a rotation amounting to 4 km 
s -1 arcsec -1 at PAs of 62°and 102°, while Fisher et al (1995) 
find no significant rotation at PAs of 31°, 35° and 125°. 

NGC 6173 and NGC 6086 both show evidence of ro- 
tation in the core, decoupled from the kinematic proper- 
ties further out. In NGC 6173 (Figure ^(b)), a linear least- 
squares fit gives a velocity gradient of —2.56 ± 0.27 km s -1 
arcsec -1 for |r| < 7 arcsec, with no significant rotation out- 
side this, and in NGC 6086 (Figure m(c)), similar fits give a 
velocity gradient of 0.86 ± 0.48 km s -1 arcsec -1 for \r\ < 6 
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NGC 6166 major axis 
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Figure 2. (a) The t>, a, /13, /14, /is & he along the major axes for NGC 6166. The la errorbars are also plotted. A weighted least-squares 
fit to the velocity profile is shown, with a slope of 1.07 km s _1 arcsec -1 . 
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NGC 6173 major axis 
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NGC 6086 major axis 
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arcsec and —1.33 ± 0.13 km s _1 arcsec -1 for 6 < \r\ < 32 
arcsec. Although the core rotation is only marginally sig- 
nificant, it does appear to be counter-rotating. NGC 6173 
is a shell galaxy, and both the shells and higher core rota- 
tion could be a result of a merger or interaction. However, 
in both galaxies, there is no significant asymmetry in the 
line profiles (i.e. changes in /13) associated with the peculiar 
kinematics, as expected from classical KDCs like IC 1459 
or NGC 2865 (Hau, Carter, & Balcells 1998). Also, Statler 
(1991) proposes that KDCs might also arise from streaming 
in a triaxial potential, without a merger/interaction. 

3.3 hi 

The parameter hi measures symmetric deviations from a 
Gaussian profile and is related to the velocity anisotropy. hi 
is significantly positive, with no radial dependence, in all of 
our galaxies. Although it is not straightforward to translate 
hi to the velocity anisotropy parameter f3 (= 1 — Og/cr^), 
positive values, at least outside the nucleus, generally indi- 
cate a bias towards radial orbits (Gerhard 1993; Gerhard et 
al. 1998, Rix et al. 1997). The most important result from 
the line-profile analysis is that there is no change in hi as- 
sociated with the rising dispersion in NGC 6166, indicating 
that the increase in velocity dispersion is not associated with 
a change in velocity anisotropy towards tangential orbits. 

h 4 has a magitude of 0.10, 0.05 & 0.05 in NGC 6166, 
6173 & 6086 respectively. Template mismatch can introduce 
systematics in hi by either increasing or decreasing its mag- 
nitude; see, for example, van der Marel et al. (1994), and 
Hau, Carter, & Balcells (1998). The x 2 fits in NGC 6173 
& NGC 6086 appear to be near perfect. Inspection of the 
residuals for NGC 6166 shows that the fits are less perfect 
than for the other two galaxies, but are satisfactory. We 
feel that template mismatch is not a major problem for all 
three galaxies as the optimal template is estimated from 28 
template stars. 

In order to estimate the effect of template mismatch on 
the absolute value of hi (and other kinematic parameters), 
in Appendix A we derive the kinematic parameters for NGC 
6166 for a range of individual templates. None of these tem- 
plates fit the galaxy data as well as the optimal template. 
Given an extreme range of templates, we find zero point dif- 
ferences of up to 0.1 in hi, /15 and h^; 0.2 in /13, and 50 
km/s in velocity dispersion. However these differences are 
largely in the zero points for stars at the extreme end of 
the range of spectral types of the observed templates (G6 
III and Ml III). In particular low hi is found only when 
an Ml III template is used, this template does not fit the 
data well and contributes little to the optimal composite 
templates. Although template mismatch may in principle 
shift the zero-point in hi, it cannot explain the significantly 
positive hi seen in our data. 

In Appendix A we also analyse the effect of noise on the 
recovery of hi when the magnitude of the real value of hi 
is large. We find that there are no systematic biases unless 
\hi\ > 0.1, when \hi\ is systematically underestimated at 
low signal to noise. Thus the only possible bias in our results 
would be to underestimate hi in the outer regions of NGC 
6166, although there is no evidence that this has happened, 
we cannot absolutely rule it out. 

Van der Marel & Franx (1993) give formulae which re- 



late hi to /3, but they are model-dependent. To better quan- 
tify the anisotropy parameter (5, we need to solve for the 
M/L ratio using all of the available kinematic and photo- 
metric data. Such techniques have been developed by Rix et 
al. (1997), using Schwarzchild's method of populating orbits, 
and by Dehnen (1995) for the axisymmetric case in order to 
model M32. 

In summary, the rising velocity dispersion profile in 
NGC 6166 is not due to a change in velocity anisotropy, 
but rather, reflects an increasing M/L ratio. This indicates 
that NGC 6166 possesses a massive dark halo. 



4 GAS KINEMATICS IN NGC 6166 

Significant [OIII] and [NI] emission is observed in the centre 
of NGC 6166. In this section we present the line-of-sight 
velocity distribution of the ionized gas inferred by the profile 
of the [OIII] and [NI] emission lines, obtained directly from 
the residual map of the Gauss-Hermite analysis, described 
in the previous section. This method is very successful and 
is more objective than fitting the baseline by hand. The 
results are plotted in Figure ^[ with the stellar rotation curve 
from Figure ^(a) over-plotted. We see that the [OIII] and 
[NI] emission is very concentrated to the galaxy centre, with 
spatial FWHM of 1.9 and 3.0 arcsec respectively. The spatial 
FWHM of template stars is about 2.3 arcsec, thus the [OIII] 
emission region is unresolved but that of the [NI] may be 
marginally resolved. At the galaxy centre, the [OIII] and 
[NI] have velocity dispersions a of 333 ± 11 and 365 ±20 km 
s _1 respectively, which is very similar to the central stellar 
velocity dispersion (Figure 2a). 



5 STELLAR ABSORPTION LINE 
ABUNDANCES 

5.1 Introduction 

In this section the absorption line indices are presented for 
the three galaxies. Mg2, Fes27o, and Fe5335 are extracted 
from individual, de-redshifted spectra following the defini- 
tions and recipes of Faber et al. (1985), and adopting the 
passband definitions of Worthey et al. (1994); see Appendix 
B for more details. We first discuss the local Mg2— <Fe> re- 
lations for each galaxy, then present the Mg2 radial profiles. 
We discuss the significance of these results in our Conclu- 
sions. 

As discussed in Appendix B, there is some uncertainty 
regarding the zero-point of our measured indices in the 
Lick/IDS scale, due to the fact that only one star with Lick 
indices (HR 5227) has been observed, and also that the grat- 
ing settings for the galaxy and the template star observa- 
tions were different. However recent observations by Trager 
et al. (1998) for the nucleii of NGC 6166 and NGC 6086 
do agree within their admittedly rather large error bars. We 
also find that our Mg2 measurements agree well with those 
measured by Cardiel et al. (1997), if both are not zero-point 
corrected (Figure |^). As most of our discussion is not based 
on the absolute values of the indices, throughout this Section 
we present line indices without zero-point corrections. 
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Figure 3. The gas kinematics from the [OIII] and [NI] emission lines at 5007A and 5200A along the major axis of NGC 6166 are ploted 
in the left and right plots respectively. The stellar velocity curves from Fig. 0(a) are over-plotted in white. 



5.2 Local <Fe> Versus Mg 2 

In Figure [| we plot <Fe> [= (Fe 52 70 + Fe 5 33 5 )/2)] against 
Mg2 for the three galaxies. All three galaxies appear to 
have [Mg/Fe] ratios higher than the model values with so- 
lar [Mg/Fe] abundance, but there are uncertainties in our 
zero-points. Notice further from Figure ^ that NGC 6086 
and NGC 6166 are more 'enhanced' than NGC 6173. This 
is interesting, given that the first two galaxies are centrally- 
located cDs, while NGC 6173 is best described as a gE in 
a subcluster in Abell 2197. Similarly, M87 is offset from the 
non-BCGs in the data of Davies, Sadler, & Peletier. 

Hau (1998) has pointed out that in galaxies of high ve- 
locity dispersion, a positive value of hi, such as might result 
from radial anisotropy, leads to more absorption being lost 
from the passband that the Fe indices are measured from, 
and thus to an artificially high value of [Mg/Fe]. This sys- 
tematic error might account for as much as a third of the off- 
set between NGC 6166 and NGC 6173. However we still find 
that our galaxies with cD morphology have higher [Mg/Fe] 
compared with NGC 6173, which is the only normal bright 
elliptical in our sample. If further observations confirm that 
NGC 6173 is typical (despite uncertainty in the relative zero- 
points, the data of Davies et al. (1993) support this), then 
there is a difference between cDs and normal bright ellip- 
ticals. Bender (1996) uses the offset of [Mg/Fe] to higher 
values in ellipticals when compared with spirals as an argu- 
ment against the idea that ellipticals form from the mergers 
of objects similar to present-day spirals; any such argument 
is clearly stronger for cDs. Bender (1996) and Davies (1996) 
feel that the most plausible explanation for the non-solar 
[Mg/Fe] ratios in ellipticals is a short burst of star formation, 
< 1—2 Gyr. This would produce a population enriched in 
Mg from rapid Type II supernovae from massive stars, while 
few stars would be formed from Fe-enriched gas formed in 
longer-lived Type la supernovae. See Section 6 for further 
discussion on these topics. 



5.3 Mg2 Abundance Gradients 

We first compare our Mg2 profile for NGC 6166 with that of 
Cardiel et al. (1997) in Figure |H| We have only plotted the 
SW side of the major axis because of contamination from 
the secondary nucleus on the NE side. Neither set of data 
has had any zeropoint corrections applied (see Appendix B). 
Out to ~ 15 arcsec, the agreement is very good. Beyond 15 
arcsec, our Mg2 is significantly higher than that of Cardiel 
et al. The source of this discrepancy is unclear, however our 
signal-to-noise is higher than that of Cardiel et al. Their 
three outermost points have Mgg values < 0.2, lower than 
any of the ellipticals studied by Davies et al. (1993). These 
extreme values seem unlikely, and we are confident of the 
reliability of our Mg2 profiles. 

The sensitivity of our line strength measurements to sky 
subtraction errors needs to be considered. The sky spectrum 
that we subtract is averaged over many spatial pixels and is 
relatively low noise. The only feature in the sky spectrum 
are the night sky emission line at 5577A and (much weaker) 
5199A, apart from this the sky is well represented by a con- 
stant continuum. At these redshifts the [NI] line at 5199A 
does not fall in the Mg2 band or its continuum bands, and so 
can be neglected. Therefore the sky in the region of all of the 
line and continuum bandpasses can be represented by a con- 
stant, and the error in this constant can be determined from 
the maximum residuals seen in the [OI] line at 5577A in the 
sky subtracted spectra, and the ratio of the peak intensity 
in this line to the continuum in the sky spectra. 

Inspection of the sky line residuals in the outer regions 
of our profile show that sky subtraction errors do not exceed 
1% of the signal level, thus systematic errors in our Mg2 
values will not exceed 0.01 mag. 

In Figure we present the major axis profiles of Mg2, 
for NGC 6166, 6173, and 6086. Outside 3-5 arcsec, the Mg 2 
profiles are well-fit by power-laws of slope (d(Mg2)/d(log 
r)) -0.080, -0.059, and -0.073 for NGC 6166, 6086 and 
6173 respectively (inside 3 arcsec, seeing effects are likely to 
be important; the flattening of Mg2 in NGC 6166 within a 
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Figure 4. Local <Fe> Versus Mg2 for three BCGs. Open triangles are plotted for NGC 6166, filled squares for NGC 6086, and open 
circles for NGC 6173. No zcropoint corrections have been applied to the data in Figures 4, 5 and 6. For NGC 6173 and NGC 6086, the 
data have been folded about the major axis, while for NGC 6166 we only show the SW side of the major axis, to avoid contamination 
from the secondary nucleus on the NE side. The dark solid lines are weighted linear least-squares fits to the data points; from left to right, 
they correspond to NGC 6173, NGC 6086 and NGC 6166. Plotted also are the locations expected for a 12 Gyr old population with 0.4 
< Z/Zq < 2.5, predicted by 3 population synthesis models: Worthey (1994; thin solid line); Bruzual & Chariot (1997; dot-dashed line), 
and Vazdekis et al. (1996; dashed line). The line fitted to the nuclear indices of a sample of galaxies from Worthey, Faber, & Gonzalez 
(1992) is indicated by the dotted line. 



larger radius of ~ 5 arcsec may be due to star formation 
inside the emission- line region in this galaxy). These slopes 
are fairly similar to each other, and are also similar to the 
average Mg2 slope of —0.058 found by Cardiel et al. (1997) 
for eight BCGs without cooling flows or optical emission 
lines. Gorgas, Efstathiou, & Aragon-Salamanca (1990) also 
measured shallower gradients for two BCGs (—0.051 and 
-0.021 for the cDs in 0559-40 and PKS 2354-35). However, 
examination of their Figure 9 shows that these gradients are 
based on only three points, and the outermost point in each 
case has large uncertainty. 

Cardiel et al. (1997) find that the Mg2 slopes are 
shallower in BCGs with optical line-emission, inside the 
emission-line region; our data for NGC 6166 are consistent 
with this. Cardiel et al. attribute this flattening to star for- 
mation from cooling flow gas near the galaxy centre. This 
makes comparison of BCG Mg2 data difficult, unless one is 
careful to compare slopes outside emission-line regions in all 
cases. 



We can also compare our Mg2 slopes with those of non- 
BCG ellipticals. Couture & Hardy (1988) found a mean 
Mg2 gradient of —0.053 ± 0.015 for six early-type galax- 
ies; Davies, Sadler, & Peletier (1993) find a mean slope of 
-0.059 ± 0.022 for 13 normal ellipticals; Gorgas et al. (1990) 
obtain a mean slope of —0.058 ± 0.027 for 16 ellipticals; 
Davidge (1992) found a steeper gradient of —0.081 ± 0.01 
for 11 ellipticals. Gonzalez & Gorgas (1998; see also Gonza- 
lez & Gorgas 1997) find a mean Mg2 gradient of —0.055 ± 
0.025 for 109 early-type galaxies. Thus, our three BCGs have 
slopes consistent with those found in non-BCG ellipticals. It 
seems likely that these BCGs have not experienced signifi- 
cantly more mergers than normal ellipticals, as mergers are 
expected to dilute abundance gradients (White 1980). 



6 CONCLUSIONS 

Our main results are: 
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Figure 5. Mg2 profiles for NGC 6166. Mg2 is plotted against radius along the SW side of the major axis. Closed circles are plotted for 
our data, while open triangles are used for the data of Cardiel et ai. (1997). Both datasets have not been zero-point corrected. 



(1) : We find small rotation along the major axes of NGC 
6166, NGC 6173, and N6086, with V rot < 50 km s" 1 at 
the outermost observed major axis radii (30-40 arcsec) in 
the three galaxies; the corresponding V/cr ranges from 0.02- 
0.18. This small rotation is consistent with the nearly com- 
plete lack of rotation found near the centres of a sample 
of 13 BCGs by Fisher et aJ. (1995). This is consistent with 
merger models of BCG formation, but data at larger radii 
are urgently needed, since numerical simulations predict 
that mergers should contain significant amounts of angular 
momentum beyond 2—3 r e (e.g. Hernquist 1993). 

NGC 6173 and NGC 6086 have larger velocity gradients 
near the galaxy centres, and NGC 6086 appears to have a 
counter-rotating core. These cores may arise from a previous 
merger/interaction (supported by the shells seen in NGC 
6173), or may be due to streaming in a triaxial potential; 
further observations at other position angles will be needed 
to decide between these two alternatives (Statler 1991). 

(2) : NGC 6086 and NGC 6173 have slowly declining ve- 
locity dispersion profiles. In contrast, NGC 6166 has a rising 
velocity dispersion: a increases from ~ 325 km s _1 at the 
centre to ~ 440 km s~ x at 35arcsec. 

(3) : The Gauss- Hermite /u moment is positive and roughly 
constant at 0.05 at all radii along the major axis of NGC 
6173 and NGC 6086, and 0.1 along the major axis of NGC 
6166. Template mismatch is unlikely to account for all of 
this offset, and we conclude that there is possibly radial 



anisotropy present. We will confirm and quantify this sug- 
gestion with detailed modelling in a later paper. 

(4) : We find [OIII] and [NI] emission concentrated within 
5 arcseconds (3h _1 kpc) of the centre of NGC 6166, with a 
central velocity dispersion of ~ 350 km s _1 , similar to the 
central stellar velocity dispersion. The extent of the [OIII] 
emission region is unresolved, but that of the [NI] may be 
marginally resolved. 

(5) : We find that [Mg/Fe] is larger in NGC 6166 and NGC 
6086 than in NGC 6173. Although positive Ha can lead to 
systematic errors in [Mg/Fe], the magnitude of the offset 
appears too large to be explained by such an error. The 
difference in [Mg/Fe] between galaxies with cD morphology 
and normal bright ellipticals, if confirmed by further obser- 
vations, would present a difficulty for models in which the 
cDs were formed from mergers of normal galaxies of any 
type. 

(6) : The Mg2 gradients outside 3-5 arcsec range from 
—0.06 to —0.08 for the three galaxies. These slopes are con- 
sistent with those found for BCGs without cooling flows or 
emission lines by Cardiel et aJ. (1997), and for non-BCG 
ellipticals. More data are needed to definitively compare 
BCGs with non-BCG ellipticals in this regard. 

We defer to a forthcoming paper a detailed dynamical 
modelling of these three galaxies, but we note that (2) and 
(3) together imply that the M/L ratio increases outwards 
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Figure 6. Mg 2 profiles for NGC 6166, NGC 6173, and NGC 6086. Mg 2 is plotted against the log of major axis radius. For NGC 6173 
and NGC 6086, data have been folded about the major axis, while for NGC 6166 only the SW side of the major axis is shown because 
of problems with the secondary nucleus on the NE side. The solid line in each plot is a weighted least-squares fit, with |r| < 3 arcsec 
excluded for NGC 6173 and NGC 6086, and \r\ < 5 arcsec excluded for NGC 6166. 
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in the NGC 6166 halo, and likely in the other two galaxies 
as well. Further support for a DM halo in NGC 6166 comes 
from X-ray data, e.g. Buote & Canizares (1996). 

Do these results favour one BCG formation model over 
others? There are five main scenarios for BCG formation: 
(i) Deposition of cooling flow gas from the ICM (e.g. Fabian 
1994); (ii) Cannibalism of other cluster galaxies (Ostriker 
& Tremaine 1975); (iii) Tidal stripping of other cluster 
galaxies, either by the cluster tidal field or by the BCG 
itself (Richstone 1976; Moore et al. 1996); (iv) Earlier 
origin in smaller groups/subclusters, which later merged to 
form the present-day cluster (Merritt 1984, 1985; Tremaine 
1990); and (v) Dissipative collapse models (Larson 1975; 
Carlberg 1984a, b). These scenarios are not necessarily mu- 
tually exclusive; for instance, many authors have suggested 
that mergers/cannibalism create the central body of the cD, 
while tidal stripping accounts for their extended halos. 

Dubinski (1998) has run an impressive simulation of 
cD formation, using a large cosmological N-body simulation 
of cluster collapse with hierarchical merging. His 'cluster' 
has a mass of 1.0 x 10 14 Mq and a line-of-sight velocity 
dispersion of 550 km s _1 , and he notes that it would be 
classified as a poor cluster or large group by observers. He 
finds that the final central giant elliptical has a small ro- 
tation of ~ 50 km s" 1 , is radially anisotropic in the outer 
regions, is aligned in position angle with the 'cluster' galaxy 
distribution, and has a slowly declining major-axis veloc- 
ity dispersion profile. These properties agree with those ob- 
served in many BCGs, including our sample described in 
this paper. All three of our galaxies have little or no rota- 
tion, and have significantly positive /14 values which gener- 
ally indicate radial anisotropy. NGC 6086 and NGC 6173 
have velocity dispersion profiles similar to Dubinski's rem- 
nant, and NGC 6173 and NGC 6166 are aligned (to within 
20 degrees) of their host cluster/subcluster galaxy distribu- 
tions and X-ray isophotes (Carter & Metcalfe 1980; Dixon, 
Godwin, & Peach 1989; Muriel, Bohringer, & Voges 1996; 
Buote & Canizares 1996). However, more such simulations 
are needed to see if, among other things, remnants with in- 
creasing dispersion profiles and extended halos can be cre- 
ated. Garijo, Athanassoula, & Garcia-Gomez (1997) have 
run a larger set of simulations of cD formation, with a va- 
riety of initial conditions. Some of the simulations produce 
objects with extended halos and radial anisotropy, and all of 
the final galaxies have declining velocity dispersion profiles. 
Note, however, that these simulations were run with all of 
the mass initially bound to galaxies, probably not physically 
realistic. The radial anisotropy found by both authors is a 
result of the 'cosmological' mergers that are being modelled, 
where the initial filamentary structure means that merging 
galaxies enter the central object on mainly radial orbits. 

Further constraints can be placed on formation models 
if one considers the stellar abundance gradients. Our three 
galaxies all show modest gradients, comparable with nor- 
mal bright ellipticals (e.g. Davies, Sadler, & Peletier 1993). 
These gradients are not as steep as predicted by simple dis- 
sipational collapse models (e.g. Larson 1975), however they 
are steeper than we would expect from hierarchical dissi- 
pationless merger models, which predict no gradients. It is 
not clear what the net effect of gaseous mergers will be, 
since on the one hand mergers are expected to dilute exist- 



ing abundance gradients (White 1980), while on the other 
hand central star formation induced by the merger could 
augment any gradient (see for instance Mihos & Hernquist 
1994). In any case the gradients point to a formation history 
not too different from that of normal bright ellipticals. 

The [Mg/Fe] ratio in NGC 6166 and NGC 6086, and 
the difference between these two BCGs and NGC 6173, pro- 
vide important new clues as to the differences between cDs 
and normal bright ellipticals. If this difference is confirmed 
by further observations, this is a second property, alongside 
the globular cluster specific frequency, in which ellipticals 
differ from spirals and in which cDs are even more extreme. 
Fisher et al. (1995) and Bender (1996) find that [Mg/Fe] 
is correlated with velocity dispersion, so perhaps the large 
[Mg/Fe] offsets seen in our BCGs are simply the contin- 
uation of this correlation to the largest measured velocity 
dispersions. The trend of increasing [Mg/Fe] we see in our 
three BCGs (Figure 4) does agree with their relative velocity 
dispersions. A dissipative process is most likely to give rise 
to these differences, although mergers as a trigger for such a 
process cannot be ruled out; for instance Zepf and Ashman 
(1993) interpret the bimodal colour (and hence metallicity) 
distributions for globular clusters in some elliptical galaxies, 
including the BCGs M87 and NGC 1399, as evidence for a 
second, merger induced, epoch of cluster formation. 

If we had data for NGC 6166 and NGC 6173 alone, it 
would be tempting to draw a clear distinction between the 
two galaxies. NGC 6166 has a rising velocity dispersion, and 
enhanced [Mg/Fe] , a high globular cluster specific frequency 
(Bridges et al. 1996) and an X-ray cooling flow. NGC 6173 
has none of these things ^, yet is a galaxy with a similar lu- 
minosity and surface brightness profile (Carter 1977). NGC 
6086 appears to have a mixture of properties: it has a high 
[Mg/Fe], yet does not have a rising dispersion profile. It 
is important to establish which of these properties depend 
upon the environment of a galaxy, and similar observations 
of galaxies with cD morphology in low density environments 
(such as NGC 4839; Oemler 1976) are required. 

It is not clear how to interpret the rising velocity dis- 
persion profile that we observe in NGC 6166. One would 
expect a rising dispersion profile under any model of BCG 
formation, since under simple energy arguments the BCG 
dispersion profile should eventually join that of the clus- 
ter as a whole (as is the case for NGC 1399/Fornax; see 
Section 1). However, we note that the three BCGs with in- 
creasing dispersions, and where X-ray data for the cluster 
exist (NGC 6166, IC 1101, Abell 3266) are all in rich clus- 
ters with smooth X-ray morphology. There is no evidence 
for substructure or peculiar BCG velocities in any of these 
clusters. We suggest that these are old, relaxed clusters in 
which the BCG is at the centre of the cluster potential; these 
are the sorts of clusters in which it will be most profitable 
to search for BCGs with increasing velocity dispersions, in 
order to study the transition between the BCG and the clus- 
ter. 



* The evidence that NGC 6173 does not have a high globular 
cluster specific frequency comes from a deep, unpublished, R band 
image taken at the same time as the R-band image of Bridges et 
al. (1996) 
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APPENDIX A: MONTECARLO SIMULATIONS 
OF THE SENSITIVITY OF THE RECOVERED 
KINEMATIC PARAMETERS TO NOISE 

In this appendix the sensitivity of the kinematic parameters 
to noise for NGC 6166, NGC 6173 & NGC 6086 is investi- 
gated via Montecarlo simulations of the effect of both Pois- 
son noise and sky residuals. Especially of interest is whether 
the lowered velocity dispersions measured for points beyond 
40arcsec for NGC 6166 and beyond 20arcsec for NGC 6173 
are significant. For this purpose artificial spectra are gen- 
erated by convolving the optimal stellar templates at the 
centres of NGC 6166 & NGC 6173 with perfect Gaussian 
LOSVDs. 

Poisson noise with different random seeds is added to 
the artificial spectra. Imperfect sky subtraction is investi- 
gated by adding blank spectra extracted at large distances 
from the nucleus. Up to 8 rows are added to mimic the worst 
case of sky contamination for NGC 6166, and 16 rows for 
NGC 6173. The kinematic parameters are recovered three 
times using kinematics, each with different Poisson noise 
seeds and sections of blank sky. 



NGC 6166 

For NGC 6166, LOSVDs with cz = 9255.3 kms -1 and 
a = 400, 450, & 500 kms -1 are modelled. The artificial spec- 
tra are scaled to continuum levels of 12000, 7200, 4200 & 
2000 counts per pixel, corresponding to distances of 0, 5, 8 
& ~ 15arcsec from the nucleus respectively. Beyond 15arcsec 
summation is required to achieve 2000 counts per pixel. The 
number of spectra to be summed is 2, 4, 5 & 8 for distances 
of approximately 22arcsec, 33arcsec, 37arcsec and 50arcsec 
from the nucleus respectively. The results of the experiments 



are plotted in Fig. Al, where the left, centre and right pan- 
els correspond to a — 400, 450 & 500 kms -1 respectively. 
For a = 400 kms -1 , the mean velocity v can be determined 
to within 2 kms - at distances less than 5arcsec from the 
nucleus. The scatter increases to 6 kms -1 at r = 15arcsec, 
13 kms -1 at r = 33arcsec, and 27 kms -1 at r = 50arcsec. In 
addition beyond 33arcsec there is a small systematic offset of 
~ 20 kms -1 towards higher velocity. For broader LOSVDs, 
the behavior is similar, except that the recovering accuracy 
is lower. 

If the velocity dispersion is ~ 400 kms -1 , it can be de- 
termined to within ~ 3 kms -1 at distances less than 5arcsec 
from the nucleus. Its scatter increases with distance. In ad- 
dition, as the contribution of the sky increases with decreas- 
ing S/N, the recovered value is systematically lower. The 
offset is about —12 kms -1 at r = 15arcsec, —40 kms -1 at 
r — 33arcsec, and —70 kms -1 at r = 50arcsec. The behavior 
is similar if the LOSVD is broader, except that the magni- 
tude of the systematic offset is larger. At the extreme, the 
recovered value for a Gaussian LOSVD with a — 500 kms -1 
at r — 50arcsec is 100 kms -1 lower. 

In general, the Gauss-Hermite moments hz and /14 can 
be recovered reasonably well if a = 400 kms -1 , to within 
~ 0.025 at distances up to 33arcsec. Their scatter increases 
with a. No obvious systematics are introduced in /13 even at 
the outermost points, or when a is increased to 500 kms -1 . 
A small positive offset of order 0.03 is introduced in hi for 
data points beyond 33arcsec, and this tends to increase with 
a. The Gauss-Hermite moments /15 and ha can be recovered 
well to within 0.01 in the inner 15arcsec, and to ~ 0.02 at 
r = 33arcsec. Again, the scatter increases with a. In addi- 
tion a significantly positive offset of ~ +0.05 is introduced 
to the hs at the outer points if sigma is increased beyond 
400 kms -1 . 

The conclusions from this exercise are that at dis- 
tances beyond ~ 33arcsec, where summation of more than 
4 spectra are required, systematics start to creep in as 
the sky contribution increases. At such distances, the re- 
covered velocity dispersions are systematically lower. This 
shows that we are probably unable to measure LOSVDs 
broader than ~ 400 kms -1 beyond these distances, and the 
low (~ 320kms -1 ) velocity dispersions beyond 40arcsec ac- 
tually measured for NGC 6166 are entirely consistent with 
expectation. For these reasons, we measure v up to 50arcsec, 
and a, /13, /14, h5 & h6 up to 36arcsec from the nucleus. 



NGC 6173 

For NGC 6173, LOSVDs with cz = 8815.3 kms -1 and a = 
200, 250, & 300 kms -1 are modelled. The artificial spectra 
are scaled to continuum levels of 17000, 3600 & 2000 counts 
per pixel, corresponding to distances of 0, 6, & ~ 9arcsec 
from the nucleus respectively. Beyond 9arcsec summation is 
required to achieve 2000 counts per pixel. The number of 
spectra to be summed is 4, 8 & 16 for distances of approx- 
imately 19arcsec, 30arcsec and 46arcsec from the nucleus 
resp ecti vely. The results of the experiments are plotted in 
where the left, centre and right panels correspond to 
250 & 300 kms -1 respectively. For a = 200 kms -1 , 
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Figure Al. The effect of noise on the recovery of the kinematic parameters for NGC 6166. The composite stellar template at the nucleus 
of NGC 6166 is convolved with different model LOSVDs to make mock galaxy spectra. Added to these are Poisson noise and blank sky 
spectra. The kinematic parameters are recovered from 3 mock spectra, each created with a different noise seed and sections of blank sky. 
The panels from left to right correspond to Gaussian LOSVDs with velocity dispersion a of 400, 450 & 500 kms -1 respectively. In each 
panel, results are shown for continuum levels of 2000, 4200, 7200 & 12000 counts per pixel, labelled along the bottom with the number 
of rows of sky summed. Up to 8 rows of sky spectra are added to the artificial spectrum to mimmick the points farthest from the galaxy 
nucleus. The model values are represented by the dotted lines, and the recovered parameter by the solid dots. 



r = 46arcsec. In addition beyond 30arcsec there is a small 
systematic offset of ~ +6kms~ 1 . For broader LOSVDs, the 
behavior is similar, except that the magnitude of the scatter 
and the systematic offset is slightly higher. 

If the velocity dispersion is ~ 250 kms" 1 , it can be de- 
termined to within ~ 6kms _1 at distances less than 6arcsec 
from the nucleus. Its scatter increases with distance. In ad- 
dition, as the contribution of the sky increases with decreas- 
ing S/N, the recovered value is systematically lower. The 
offset is about — 9kms _1 at r = 19arcsec, — 15kms _1 at 
r — 30arcsec, and — 27kms _1 at r — 46arcsec. The mag- 



nitude of the systematic offset is larger if the LOSVD is 
broader, and smaller if the LOSVD is narrower. At the 
extreme, the recovered value for a Gaussian LOSVD with 
a = 300 kms -1 at r — 46arcsec is 48kms~ 1 lower. 

In general, the Gauss-Hermite moments and /14 
can be recovered reasonably well if the S/N is high. If 
a — 250 kms -1 , they can be determined to within ~ 0.02 at 
distances up to 9arcsec. Their scatter increases with a and 
sky contribution. A small positive offset in /13 is introduced 
with increasing sky contribution. No significant systemat- 
ics are introduced in /14, except for the outermost points 
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Figure A2. Same as Fig. |Al| , but for NGC 6173. The panels from left to right correspond to Gaussian LOSVDs with velocity dispersion 
a of 200, 250 & 300 kms -1 respectively. In each panel, results are shown for continuum levels of 2000, 3600 & 17000 counts per pixel, 
labelled along the bottom with the number of rows of sky summed. Up to 16 rows of sky spectra are added to the artificial spectrum to 
mimmick the points farthest from the galaxy nucleus. The model values are represented by the dotted lines, and the recovered parameter 
by the solid dots. 



ft 



4 • 



at r ~ 46arcsec, where hi tends to be lower. The Gauss- 
Hermite moments /15 and h§ can be recovered well only 
near the galaxy centre (r < 9arcsec). The scatter increases 
with a and increasing sky contribution. If a = 250 kms -1 , 
there are no obvious systematics in /15 with increasing sky 
contribution. However, at distances larger than 19arcsec /is 
is scattered upwards if a = 200 kms -1 , and downwards if 
a = 300 kms -1 . ha tends to be scattered downwards at dis- 
tances larger than 19arcsec for all a. 

The tests for NGC 6173 should also be applicable to 
NGC 6086, as both galaxies have similar velocity disper- 
sions, and are at similar distances. For NGC 6086, the 
continuum levels are 14600 and 2000 counts per pixel at 



r = Oarcsec and 5arcsec. Beyond 5arcsec summation is re- 
quired. The number of spectra to be summed is 4, 6 & 15 for 
distances of approximately 12arcsec, 15arcsec and 30arcsec 
from the nucleus respectively. 

The conclusions from this exercise are that at distances 
beyond ~ 19arcsec, where summation of more than 4 spec- 
tra are required, systematics start to creep in as the sky 
contribution increases. At such distances, the recovered ve- 
locity dispersions are systematically lower, but the problem 
is not as severe as in NGC 6166. This shows that the fall in a 
beyond 20arcsec is probably more gradual in reality. The sig- 
nificance of the Monte Carlo tests is that the a in NGC 6173 
is consistent with being flat or falling beyond 20arcsec, but 
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inconsistent with rising outwards. We measure v, a & ft.3 up 
to 40arcsec, hi up to 25arcsec, and /15 & /i6 up to 19arcsec 
from the nucleus. 



Recovery of High Values of hi 

The tests described above show that the systematic errors 
on the recovery of hi caused by noise are small in the case 
of /14 = 0. It is important also to ask what other values of hi 
can be recovered by our techniques. For NGC 6166 we model 
LOSVDs with cz = 9255.3 kms" 1 ans a = 400 kms" 1 , and 
hi =-0.10, -0.05, 0.00, +0.05, +0.10 and +0.15 respectively, 
and the results are plotted in Fig. A3. The systematic bias 
on a at low S/N is marginally worse for positive hi (second 
panel up in each column), and the main result of this exer- 
cise is that there is a systematic underestimation of |/u| for 
\hi\ > 0.1. Our measurement of high hi in NGC 6166 could 
possibly still be a systematic underestimate. 



The Effect of Template Choice 

Our kinematic parameters are all determined using optimal 
best fit templates constructed from 28 stellar spectra. To 
demonstrate that this procedure does not itself introduce 
systematic effects into the kinematic parameters, here we 
repeat the analysis for NGC 6166 using individual template 
stars. Moreover to demonstrate the effect of template choice 
on the derived parameters we have carried out this exercise 
with templates with a wide range of spectral type. In Fig. A.4 
we plot the kinematic parameters for NGC 6166 derived 
from the five individual templates: HR7555 (G6III, open 
squares); HR4864 (G7V, filled triangles); HR6674 (K0III, 
crosses); HR4676 (K4III, open circles) and HR8057 (M1III, 
thin rectangles). The trends in the derived parameters are 
identical to those in Fig. |^, although there are differences 
in the zero points which depend upon the template choice. 
The difference in the velocity zero points is not real, it is 
because this procedure determines the velocity relative to 
the first star in the list of templates, if there is only one it 
depends upon the radial velocity of that star. For a and all 
of the Hermite terms the M1III star HR8057 gives a very 
different zero point, the dispersions are so much lower that 
man y of them are off the bottom of the dispersion panel in 
Fig. A4. However this star is a very poor fit to the galaxy 



spectra and does not contribute much to the optimal com- 
posite templates. For the other templates, hi is tightly de- 
fined around +0.10, suggesting that template mismatch er- 
rors do not contribute to our large positive values of hi. On 
the other hand there are large zero point differences in h% 
between the G and K star templates, to get reliable values 
of hz it is important that the templates used fit well. 

The velocity dispersion zero points differ with a total 
spread of 50 kms -1 (again neglecting HR8057), so the un- 
certainty of the zero point of the velocity dispersion profile 
due to template mismatch might be estimated at 25 kms -1 . 



APPENDIX B: PROCEDURE FOR 
MEASURING LICK ABSORPTION LINE 
INDICES 

In this appendix we describe the procedures for obtaining 
absorption line indices in the Lick system. Line indices are 
extracted from individual, de-redshifted spectra following 
the procedures outlined in Worthey & Ottaviani (1997), 
adopting the passband definitions of Worthey et al. (1994). 
The blue-ward continuum bandpass for Mg2 is redefined to 
be 4931.625-4957.625 A instead of 4895.125-4957.625 A, 
because the spectra of the template stars only start from 
about 4931 A. The difference between the Mg index Mg 2 ait 
under this alternative definition, and the Lick Mg2 index is 
tiny (see below). The molecular band Mg2 a!t is measured 
in pseudo magnitudes, whilst the atomic absorption line in- 
dices Fe527o and Fes335 are measured as equivalent widths 
in Angstroms. 

First, our spectra are broadened to a resolution of 
8.5A FWHM, in common with the Lick/IDS resolution in 
this part of the spectrum. In NGC 6166, there is signifi- 
cant [NI] emission at — 8.7arcsec < r < 8.1arcsec. As the 
[NI] emission line is halfway inside the Mg2 passband, its 
flux can artificially lower the Mg2 by 7% at the nucleus. 
Therefore the flux contribution between rest wavelengths 
5183.125 A and 5196.526 A is measured from the residuals 
of the kinematic fitting and subtracted off the flux in the 
Mg2 passband. 

The indices are then corrected to a zero dispersion sys- 
tem by applying correction factors. These are estimated by 
broadening the optimal stellar template of NGC 6166 to dif- 
ferent velocity dispersions, and are summarized in table B2. 
The correction for Mg2 is insignificant even at large disper- 
sions, whilst Fe5335 requires a rather large correction. These 
corrections are not sensitive to the choice of template used 
to determine these corrections, in the worst case (Fes335 at 
high velocity dispersion) the total range in the factor deter- 
mined using a range of templates is 0.034, implying an error 
in the factor as determined from the best fit template of at 
most 0.01. 



Bl Conversion from alternative system to Lick 
system 

To convert Mg2 a!t to Mg2, a scaling relationship is required. 
This is obtained by comparing the Mg2 indices for 9 tem- 
plate stars measured under the Lick and our alternative def- 
inition. The template stars were observed on the AAT in 
March 1996 using the RGO spectrograph, and their spectra 
span the wavelength range 4850-5610 A, and are broadened 
to 8.5A FWHM resolution. The fact that these spectra are 
taken with a completely different instrument and telescope 
does not pose a problem, as the only variable is the change 
in bandpass definition. The Mg2 indices under Worthey's 
definition are plotted against those under the alternative 
definition in Fig. Bl, and the following relationship is ob- 



tained by a linear least-squares fit: 



Mg 2 



1.005A/#2 



(Bl) 



Because the changes introduced by redefining the continuum 
bandpass are only 0.5% for Mg2, they are unlikely to be a 
major source of error. 



© 1998 RAS, MNRAS 000, 000-000 



Kinematics, Abundances, and Origin of BCGs 19 



0.05 
' 



0.05 
c ra 
-0.05 



-0.1 

0.05 
1 



1 400 
b 550 









T 




1 


l T 


!{ 






! t ' 


t 


t 


■ T 




i 1 t ! 

r T * 


t 


T T tj 


Mm** 


t f t t 


It If ft 




t * t * • 




| M It 


M 


t4+J 












t 










+ 

♦ 1 i . 


+ 

T + 


! M ♦ * • 




T 


tfttM 


II' 1 " 




t j t M ♦ 


H t H * 














t T 




t 


f|l ft 


j |i t 1 f 














, » • • • 

- ♦ 


♦ t • • • 


♦ J • • • 


f ♦ 1 t 

♦ 1 ♦ 

♦ 


t 1 ' « * 
* 

t 


. ! " ! 

i 


:» f 1 t , ..... 


+ 1 i * * * 


I ♦ 


* ♦ 


{ ♦ 


t ♦ 










- r -M-t t i 





^9300 h+ 
in 
\ 

I 9250 
> 

9200 r 

Counts/1 000= 2 2 2 4.2 7.2 12 2 2 2 4.2 7.2 12 2 2 2 4.2 7.2 12 2 2 2 4.2 7.2 12 2 2 2 4.2 7.2 12 2 2 2 4.2 7.2 12 
Sum= 841111841111841111841111841111841111 

Figure A3. Same as Fig. [Al] , but the six columns now represent different input values of hi, each for a = 400kms — 1 . The panels from 
left to right correspond to LO SVDs with h.4 of -0.10, -0.05, 0.00, +0.05, +0.10 and +0.15 respectively. All other parameters are as for 
the left hand panel of Fig. |Al[ 



Table Bl. The bandpass wavelengths of absorption line indices Mg2, Fe527o, and Fe5335 defined by Worthey et al. (1994), with our 
alternative definition for Mg2 aK . The latter is measured and transformed to Mg2 via the scaling relation described in the text. 



Feature Type Central Bandpass Continuum Bandpasses 

(A) (A) 

Mg 2 a " Molecular band 5154.125-5196.625 4931.625-4957.625, 5301.125-5366.125 

Mg 2 Molecular band 5154.125-5196.625 4895.125-4957.625, 5301.125-5366.125 

Fe 52 7o Atomic absorption line 5245.650-5285.650 5233.150-5248.150, 5285.650-5318.150 

Fe 53 35 Atomic absorption line 5312.125-5352.125 5304.625-5315.875, 5353.375-5363.375 



B2 Lick zero points 

Because our spectra are not flux calibrated, there might be a 
zero point difference between our measured indices and the 
Lick values, due to the curvature of the response function. Of 
our standard stars, only one, HR 5277 has been measured 
on the Lick system. Worthey (private communication) as- 
cribes low weight to the observations of this particular star, 
and in any case using observations of a standard star to de- 
termine the zero points for galaxies is subject to systematic 



error, due to the substantial redshifts of our galaxies, and 
the corresponding difference in the grating settings for the 
observations. Thus we are unable to tie our measurements 
to the Lick system directly. On the other hand Cardiel et al. 
(1997) have measurements which we, and they, believe to be 
on the Lick system. Figure |B| demonstrates good agreement 
between the Mg2 values in the region in which the data of 
Cardiel et al. are reliable, and thus we believe that the zero 
point offset between our values and the Lick system is small, 
at least in Mg2. 
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Comparison of results using different templates 
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Figure A4. Kinematic parameters for NGC 6166 derived using individual stellar template spectra: HR7555 (G6III, open squares); 
HR4864 (G7V, filled triangles); HR6674 (K0III, crosses); HR4676 (K4III, open circles) and HR8057 (M1III, thin rectangles). 



Trager et al. (1998) publish Lick indices for a 1.4 x 4 
arcsecond aperture centred on the nucleii of NGC 6086 and 
NGC 6166. For NGC 6086 they find Mg 2 = 0.309 ± 0.008; 
Fe 52 70 = 3.32 ± 0.32A; Fe 53 35 = 2.71 ± 0.47A. For NGC 
6166 they find Mg 2 = 0.315 ±0.010; Fe 52 ro = 3.21±0.34A; 
^65335 = 2.01 ± 0.47 A These are consistent with our re- 
sults, although the error bars on their iron line equivalent 
widths are disappointingly large. Although some uncertainty 
remains as to the zero point of the iron line indices, it seems 
that our indices are within 0.01 magnitudes of the Lick sys- 
tem in the case of Mg2, and within 0.3 Angstrom in the case 
of the combined iron line index. 
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Kinematics, Abundances, and Origin of BCGs 



Table B2. Velocity dispersion correction factors for different line 
indices. The measured indices should be divided by these to trans- 
form back to the zero dispersion system. 



cr Mg 2 
(kms -1 ) 

160 0.998 

180 0.998 

200 0.997 

220 0.997 

240 0.997 

260 0.996 

280 0.996 

300 0.995 

320 0.995 

340 0.994 

360 0.993 

380 0.991 

400 0.990 

420 0.988 

440 0.987 



Fe 52 7o Fe 53 35 



0.938 0.895 

0.924 0.868 

0.909 0.842 

0.895 0.815 

0.879 0.785 

0.864 0.756 

0.849 0.727 

0.834 0.697 

0.819 0.668 

0.804 0.640 

0.788 0.612 

0.773 0.586 

0.758 0.560 

0.742 0.536 

0.726 0.512 



Mg2 



a 

— • o 
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Figure Bl. Mg2 measured according to the Lick definition is 
plotted against that measured according to our alternative defi- 
nition. 
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